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ABSTRACT 

These lectures provide an introduction to lattice methods for nonperturbative stud- 
ies of Quantum Chromodynamics. Lecture 1 (Ch. 2): Basic techniques for QCD 
and results for hadron spectroscopy using the simplest discretizations; lecture 2 (Ch. 
3): "improved actions" -what they are and how well they work; lecture 3 (Ch. 4): 
SLAC physics from the lattice: structure functions, the mass of the glueball, heavy 
quarks and a s (Mz), and B — B mixing. 



1 Introduction 

The latticed version of QCD was invented by Wilsoni in 1974. It has been a fruitful source 
of qualitative and quantitative information about QCD, the latter especially in the years since 
Creutz, Jacobs, and Rebbii performed the first numerical simulations of a lattice gauge theory. 
Lattice methods are presently the only way to compute masses and matrix elements in the strong 
interactions beginning with the Lagrangian of QCD and including no additional parameters. In 
the past few years the quality of many lattice predictions has become very high, and they are 
beginning to have a large impact in the wide arena of "testing the standard model." My goal 
in these lectures is to give enough of an overview of the subject that an outsider will be able to 
make an intelligent appraisal of a lattice calculation when s/he encounters one later on. 

The first lecture will describe why one puts QCD on a lattice, and how it is done. This is a 
long story with a lot of parts, but at the end I will show you "standard" lattice results for light 
hadron spectroscopy. The main problem with these calculations is that they are so unwieldy: to 
get continuum-like numbers requires very large scale numerical simulations on supercomputers, 
which can take years to complete (sort of like the high energy experiments themselves, except 
that we do not have to stack lead bricks). We would like to reduce the computation burden of 
our calculations. In Lecture Two I will describe some of the different philosophies and techniques 
which are currently being used to invent "improved actions." Some of these methods actually 
work: some QCD problems can be studied on very large work stations. Finally, in Lecture Three 

1 will give a survey of recent lattice results for matrix elements, using physics done at SLAC as 
my unifying theme. 

2 Gauge Field Basics 
2.1 Beginnings 

The lattice is a cutoff which regularizes the ultraviolet divergences of quantum field theories. As 
with any regulator, it must be removed after renormalization. Contact with experiment only 
exists in the continuum limit, when the lattice spacing is taken to zero. 

We are drawn to lattice methods by our desire to study nonperturbative phenomena. Older 
regularization schemes are tied closely to perturbative expansions: one calculates a process to 
some order in a coupling constant; divergences are removed order by order in perturbation theory. 
The lattice, however, is a nonperturbative cutoff. Before a calculation begins, all wavelengths less 
than a lattice spacing are removed. Generally one cannot carry out analytical studies of a field 
theory for physically interesting parameter values. However, lattice techniques lend themselves 
naturally to implementation on digital computers, and one can perform more-or-less realistic 
simulations of quantum field theories, revealing their nonperturbative structure, on present day 
computers. I think it is fair to say that little of the quantitative results about QCD which have 
been obtained in the last decade, could have been gotten without the use of numerical methods. 

On the lattice we sacrifice Lorentz invariance but preserve all internal symmetries, including 
local gauge invariance. This preservation is important for nonperturbative physics. For example, 
gauge invariance is a property of the continuum theory which is nonperturbative, so maintaining 



it as we pass to the lattice means that all of its consequences (including current conservation 
and renormalizability) will be preserved. 

It is very easy to write down an action for scalar fields regulated by a lattice. One just 
replaces the space-time coordinate by a set of integers (x^ — an M , where a is the lattice 
spacing). Field variables 4>(x) are defined on sites 4>{x n ) = <j) n , The action, an integral over the 
Lagrangian, is replaced by a sum over sites 
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and the generating functional for Euclidean Green's functions is replaced by an ordinary integral 
over the lattice fields 

Z= [([[^e-P 3 . (2) 

J n 

Gauge fields are a little more complicated. They carry a space-time index [i in addition to an 
internal symmetry index a [A a ^(x)) and are associated with a path in space x M (s): a particle 
traversing a contour in space picks up a phase factor 

ip — > P(exp ig J dx^A^ip (3) 

= U{s)^{x). (4) 

P is a path-ordering factor analogous to the time-ordering operator in ordinary quantum me- 
chanics. Under a gauge transformation g, U(s) is rotated at each end: 

U(s) -> g- 1 (x tt (s))U{s)g(x li (0)). (5) 

These considerations led Wilsoni to formulate gauge fields on a space-time lattice, as follows: 

The fundamental variables are elements of the gauge group G which live on the links of a 
four-dimensional lattice, connecting x and x + fi: U fl (x), with U^{x + fj,y = U^{x) 

t/» = exp(i fl aTM») (6) 

for SU(N). (g is the coupling, a the lattice spacing, the vector potential, and T a is a group 
generator). 

Under a gauge transformation link variables transform as 

Up(x)^V(x)U»{x)V{x + itf (7) 

and site variables as 

ip(x) -> V(x)ip(x) (8) 

so the only gauge invariant operators we can use as order parameters are matter fields connected 
by oriented "strings" of U's (Fig. la) 

ip(xi)U fl (x 1 )U fl (x 1 +fi)...ip(x2) (9) 

or closed oriented loops of U's (Fig. lb) 

Tr . . . U^U^x + /})...—> Tr .. . U^x)V^x + p)V(x + ftU^x + £).... (10) 



An action is specified by recalling that the classical Yang-Mills action involves the curl of 
A^, Fn v - Thus a lattice action ought to involve a product of U^'s around some closed contour. 
There is enormous arbitrariness at this point. We are trying to write down a bare action. So 
far, the only requirement we want to impose is gauge invariance, and that will be automatically 
satisfied for actions built of powers of traces of U's around closed loops, with arbitrary coupling 
constants. If we assume that the gauge fields are smooth, we can expand the link variables 
in a power series in gaA'^s. For almost any closed loop, the leading term in the expansion 
will be proportional to F 2 V . We might want our action to have the same normalization as the 
continuum action. This would provide one constraint among the lattice coupling constants. 

The simplest contour has a perimeter of four links. In SU(N) 

IN 

PS=— EE RcTr (l-Up{n)U v (n + fi)Ut(n + i>)Ut(n)). (11) 

This action is called the "plaquette action" or the "Wilson action" after its inventor. The lattice 
parameter (3 = 2N/ g 2 is often written instead of g 2 — 4ira s . 

Let us see how this action reduces to the standard continuum action. Specializing to the 
U(l) gauge group, and slightly redefining the coupling, 

S= ^EE Rc ( 1_ cxp(i 5 aL4 M (n) + A v {n + fi) - A„(x + v) - A v (n)})). (12) 

The naive continuum limit is taken by assuming that the lattice spacing a is small, and Taylor 
expanding 

A^n + v) =A ft (n) + ad v A ft (n) + . . . (13) 

so the action becomes 



P S = "4 E E 1 - Rc (exp(*5a[o(ft,^ - + 0(a 2 )})) (14) 

= i^ 4 EE-^ + --- as) 

d A xF 2 „ (16) 
transforming the sum on sites back to an integral. 



2.2 Relativistic Fermions on the Lattice 

Defining fermions on the lattice presents a new problem: doubling. The naive procedure of 
discretizing the continuum fermion action results in a lattice model with many more low energy 
modes than one originally anticipated. Let's illustrate this with free field theory. 
The free Euclidean fermion action in the continuum is 

S = J d 4 x$(x)^d^(x) + m^(x)^(x)]. (17) 



(a) (b) 



Figure 1: Gauge invariant observables are either (a) ordered chains ("strings") of links con- 
necting quarks and antiquarks or (b) closed loops of link variables. 

One obtains the so-called naive lattice formulation by replacing the derivatives by symmetric 
differences: we explicitly introduce the lattice spacing a in the denominator and write 

gnaive = £ - ^) + m^Mn- (18) 

The propagator is: 

G{p) = («7m sm Pfl a + ma) = 7 ^— — (19) 

2^ sin p^a + m z a z 

We identify the physical spectrum through the poles in the propagator, at p = iE: 

sinh 2 Ea = ^2 sin2 Pj a + m2fl2 ( 20 ) 
3 

The lowest energy solutions are the expected ones at p = (0, 0, 0), E ~ ±m, but there are other 
degenerate ones, at p = (n, 0, 0), (0, n, 0, ), ... (7r, 7r, 7r). This is a model for eight light fermions, 
not one. 

(a) Wilson Fermions 

There are two ways to deal with the doublers. The first way is to alter the dispersion 
relation so that it has only one low energy solution. The other solutions are forced to E ~ 1/a 
and become very heavy as a is taken to zero. The simplest version of this solution (and almost 
the only one seen in the literature until recently) is due to Wilson: add a second-derivative-like 
term 

^ = -^^^(^-2^ + ^) (21) 

to S nalve . The parameter r must lie between and 1; r = 1 is almost always used and "r = 1" 
is implied when one speaks of using "Wilson fermions." The propagator is 

_ -iln sinjyfl + ma-r £ M (cosp M Q - 1) 
2^ sin p^ a + {ma- (cos p^a - I)) 2 



It has one pair of poles at p^ ~ (±im, 0,0,0), plus other poles at p ~ r/a. In the continuum 
these states become infinitely massive and decouple (although decoupling is not trivial to prove). 

With Wilson fermions it is conventional not to use not the mass but the "hopping parameter" 
k = \{ma + 4r) _1 , and to rescale the fields tp ~ * V2n^p. The action for an interacting theory is 
then written 

S = ^™^ n ~~ K ^(^n(r - 7p) C/ A'( n )^n+A' + ^™( r + 7/i)^/£^n-|i)- (23) 
n n^i 

Wilson fermions are closest to the continuum formulation- there is a four component spinor on 
every lattice site for every color and/or flavor of quark. Constructing currents and states is just 
like in the continuum. 

However, the Wilson term explicitly breaks chiral symmetry. This has the consequence that 
the zero bare quark mass limit is not respected by interactions; the quark mass is additively 
renormalized. The value of k c , the value of the hopping parameter at which the pion mass 
vanishes, is not known a priori before beginning a simulation; it must be computed. This 
is done in a simulation involving Wilson fermions by varying k and watching the pion mass 
extrapolate quadratically to zero as ~ k c — n (n c — k is proportional to the quark mass for 
small m q ). For the lattice person, this is unpleasant since preliminary calculations are required 
to find "interesting" n values. For the outsider trying to read lattice papers, it is unpleasant 
because the graphs in the lattice paper typically list k, and not quark (or pion) mass, so the 
reader does not know "where" the simulation was done. Note also that the relation between k 
and physical mass changes with lattice coupling j3. 
(b) Staggered or Kogut-Susskind Fermions 

In this formulation one reduces the number of fermion flavors by using one component 
"staggered" fermion fields rather than four component Dirac spinors. The Dirac spinors are 
constructed by combining staggered fields on different lattice sites. Staggered fermions preserve 
an explicit chiral symmetry as m q — > even for finite lattice spacing, as long as all four fla- 
vors are degenerate. They are preferred over Wilson fermions in situations in which the chiral 
properties of the fermions dominate the dynamics-for example, in studying the chiral restora- 
tion/deconfinement transition at high temperature. They also present a computationally less 
intense situation from the point of view of numerics than Wilson fermions, for the trivial reason 
that there are less variables. However, flavor symmetry and translational symmetry are all mixed 
together. Construction of meson and baryon states (especially the A) is more complicated than 
for Wilson fermions .i 

2.3 Enter the Computer 

A "generic" Monte Carlo simulation in QCD breaks up naturally into two parts. In the "config- 
uration generation" phase one constructs an ensemble of states with the appropriate Boltzmann 
weighting: we compute observables simply by averaging TV measurements using the field vari- 
ables <ft^ appropriate to the sample 

1 N 

(r)~f = _£rfeW]. (24) 

i=l 



As the number of measurements N becomes large the quantity T will become a Gaussian dis- 
tribution about a mean value. Its standard deviation istj 



1 1 N 

*r = ^£|r[^]| 2 -r 2 ). (26) 

i=l 

The idea of essentially all simulation algorithms is that one constructs a new configuration of 
field variables from an old one. One begins with some simple field configuration and monitors 
observables while the algorithm steps along. After some number of steps, the value of observables 
will appear to become independent of the starting configuration. At that point the system is 
said to be "in equilibrium" and Eq. [24| can be used to make measurements. 

The simplest method for generating configurations is called the Metropolis! algorithm. It 
works as follows: From the old configuration {<f>} with action f3S, transform the variables (in 
some reversible way) to a new trial configuration {4>}' and compute the new action f3S' . Then, 
if S" < S make the change and update all the variables; if not, make the change with probability 
exp(-p(S> - S j). 

Why does it work? In equilibrium, the rate at which configurations i turn into configurations 
j is the same as the rate for the back reaction j — ► i. The rate of change is (number of 
configurations) x (probability of change) . Assume for the sake of the argument that Si < Sj . 
Then the rate i ^ j is NiP(i — > j) with P(i — ► j) = exp(— — Si) and the rate j — ^ i is 
NjP(j -> i) with P(j i) = 1. Thus Ni/Nj = exp(-/3(S , J - Sj)). 

If you have any interest at all in the techniques I am describing, you should write a little 
Monte Carlo program to simulate the two-dimensional Ising model. Incidentally, the favorite 
modern method for pure gauge models is overrelaxation.0 

One complication for QCD which spin models don't have is fcrmions. The fermion path 
integral is not a number and a computer can't simulate fermions directly. However, one can 
formally integrate out the fermion fields. For rif degenerate flavors of staggered fermions 

, »/ 
Z= [dU][dil>][dip]ex.p(-pS(U) -^xpMijj) (26) 

[dU](detM) n f/ 2 exp(-f3S(U)). (27) 

(One can make the determinant positive-definite by wr iting it as det(AftM) n // 4 .) The deter- 
minant introduces a nonlocal interaction among the Z7's: 

Z = J [dU] exp(-l3S(U) - ^Trln(AftM)). (28) 

All large scale dynamical fermion simulations today generate configurations using some vari- 
ation of the microcanonical ensemble. That is, they introduce momentum variables P conjugate 
to the f/'s and integrate Hamilton's equations through a simulation time t 

U = iPU (29) 

(30) 



The integration is done numerically by introducing a timestep At. The momenta are repeatedly 
refreshed by bringing them in contact with a heat bath and the method is thus called Refreshed 
or Hybrid Molecular Dynamics.! 

For special values of n/ (multiples of 2 for Wilson fermions or of 4 for staggered fcrmions) 
the equations of motion can be derived from a local Hamiltonian and in that case At systematics 
in the integration can be removed by an extra Metropolis accept/reject step. This method is 
called Hybrid Monte Carlo.i 

The reason for the use of these small timestep algorithms is that for any change in any of the 
U's, (M^M)" 1 must be recomputed. When Eq. ^ is integrated all of the U's in the lattice are 
updated simultaneously, and only one matrix inversion is needed per change of all the bosonic 
variables. 

The major computational problem dynamical fermion simulations face is inverting the fermion 
matrix M. It has eigenvalues with a very large range- from 2tt down to m q a- and in the physi- 
cally interesting limit of small m q the matrix becomes ill-conditioned. At present it is necessary 
to compute at unphysically heavy values of the quark mass and to extrapolate to m q = 0. The 
standard inversion technique today is one of the variants of the conjugate gradient algorithm.^] 

2.4 Taking The Continuum Limit, and Producing a Number in MeV 

When we define a theory on a lattice the lattice spacing a is an ultraviolet cutoff and all the 
coupling constants in the action are the bare couplings defined with respect it. When we take a 
to zero we must also specify how g(a) behaves. The proper continuum limit comes when we take 
a to zero holding physical quantities fixed, not when we take a to zero holding the couplings 
fixed. 

On the lattice, if all quark masses are set to zero, the only dimensionful parameter is the 
lattice spacing, so all masses scale like 1/a. Said differently, one computes the dimensionless 
combination am(a). One can determine the lattice spacing by fixing one mass from experiment. 
Then all other dimensionful quantities can be predicted. 

Now imagine computing some masses at several values of the lattice spacing. (Pick several 
values of the bare parameters at random and calculate masses for each set of couplings.) Our 
calculated mass ratios will depend on the lattice cutoff. The typical behavior will look like 

(ami(a))/(am 2 (a)) = mi(0)/m 2 (0) + 0(m x a) + 0((m ia ) 2 ) + ... (31) 

The leading term does not depend on the value of the UV cutoff, while the other terms do. The 
goal of a lattice calculation (like the goal of almost any calculation in quantum field theory) is 
to discover the value of some physical observable as the UV cutoff is taken to be very large, so 
the physics is in the first term. Everything else is an artifact of the calculation. We say that a 
calculation "scales" if the a— dependent terms in Eq. [H] are zero or small enough that one can 
extrapolate to a = 0, and generically refer to all the a— dependent terms as "scale violations." 

We can imagine expressing each dimensionless combination am(a) as some function of the 
bare coupling(s) {g(a)}, am — f{{g{a)}). As a — > we must tune the set of couplings {g(a)} so 

lim — f({g(a)}) — + constant. (32) 
o— >o a 



From the point of view of the lattice theory, we must tune {g} so that correlation lengths 1/ma 
diverge. This will occur only at the locations of second (or higher) order phase transitions in 
the lattice theory. 

Recall that the /^-function is defined by 

dg(a) dg(a) 
P{9) = «— = dMl/Aa)- (33) 

(There is actually one equation for each coupling constant in the set. A is a dimensional param- 
eter introduced to make the argument of the logarithm dimcnsionlcss.) At a critical point (3(g c ) 
= 0. Thus the continuum limit is the limit 

lim{. 9 (a)} {g c }. (34) 

Continuum QCD is a theory with one dimensionless coupling constant. In QCD the fixed point 
is g c = so we must tune the coupling to vanish as a goes to zero. 

Pushing this a little further, the two-loop /3-function is prescription independent, 

0(9) = -big 3 + b 2 g 5 , (35) 

and so if we think that the lattice theory is reproducing the continuum, and if we think that the 
coupling constant is small enough that the two- loop beta- function is correct, we might want to 
observe perturbative scaling, or "asymptotic scaling", m/A fixed, or a varying with g as 

aA = (^—p/^ cxp(--4 rT ). (36) 

Asymptotic scaling is not scaling. Scaling means that dimensionless ratios of physical ob- 
servables do not depend on the cutoff. Asymptotic scaling involves perturbation theory and the 
definition of coupling constants. One can have one without the other. (In fact, one can always 
define a coupling constant so that one quantity shows asymptotic scaling.) 

And this is not all. There are actually two parts to the problem of producing a number to 
compare with experiment. One must first see scaling. Then one needs to set the scale by taking 
some experimental number as input. A complication that you may not have thought of is that the 
theory we simulate on the computer is different from the real world. For example, a commonly 
used approximation is called the "quenched approximation" : one neglects virtual quark loops, 
but includes valence quarks in the calculation. The pion propagator is the propagator of a qq 
pair, appropriately coupled, moving in a background of gluons. This theory almost certainly does 
not have the same spectrum as QCD with six flavors of dynamical quarks with their appropriate 
masses. (In fact, an open question in the lattice community is, what is the accuracy of quenched 
approximation.) Using one mass to set the scale from one of these approximations to the real 
world might not give a prediction for another mass which agrees with experiment. We will see 
examples where this is important. 

2.5 Spectroscopy Calculations 

"In a valley something like a race took place. A little crowd watched bunches of cars, each 
consisting of two 'ups' and a 'down' one, starting in regular intervals and disappearing in about 



the same direction. 'It is the measurement of the proton mass,' commented Mr. Strange, 'they 
have done it for ages. A very dull job, I am glad I am not in the game.' "0 

Masses are computed in lattice simulations from the asymptotic behavior of Euclidean-time 
correlation functions. A typical (diagonal) correlator can be written as 

C(t) = (0|O(t)O(0)|0>. (37) 

Making the replacement 

0{t) = e m 0e- m (38) 
and inserting a complete set of energy eigenstates, Eq. (3.1) becomes 

CW=^|(0|O|n)| 2 e- £ " t . (39) 

n 

At large separation the correlation function is approximately 

C(t) ~ |<0jOjl>j 2 e- jBl * (40) 

where E\ is the energy of the lightest state which the operator O can create from the vacuum. 
If the operator does not couple to the vacuum, then in the limit of large t one hopes to to find 
the mass E\ by measuring the leading exponential falloff of the correlation function, and most 
lattice simulations begin with that measurement. If the operator O has poor overlap with the 
lightest state, a reliable value for the mass can be extracted only at a large time t. In some cases 
that state is the vacuum itself, in which E\ = 0. Then one looks for the next higher state-a 
signal which disappears into the constant background. This makes the actual calculation of the 
energy more difficult. 

This is the basic way hadronic masses are found in lattice gauge theory. The many calcula- 
tions differ in important specific details of choosing the operators 0(t). 



2.6 Recent Results 

Today's supercomputer QCD simulations range from 16 3 x 32 to 32 3 x 100 points and run from 
hundreds (quenched) to thousands (full QCD) of hours on the fastest supercomputers in the 
world. 

Results are presented in four common ways. Often one sees a plot of some bare parameter 
vs. another bare parameter. This is not very useful if one wants to see continuum physics, 
but it is how we always begin. Next, one can plot a dimensionless ratio as a function of the 
lattice spacing. These plots represent quantities like Eq. ||l]. Both axes can show mass ratios. 
Examples of such plots are the so-called Edinburgh plot (mjv /m p vs. m^/mp) and the Rome plot 
{mj^/m p vs. {m^/mp) 2 ). These plots can answer continuum questions (how does the nucleon 
mass change if the quark mass is changed?) or can be used to show (or hide) scaling violations. 
Plots of one quantity in MeV vs. another quantity in MeV are typically rather heavily processed 
after the data comes off the computer. 

Let's look at some examples of spectroscopy, done in the "standard way," with the plaquette 
gauge action and Wilson or staggered quarks. I will restrict the discussion to quenched simula- 
tions because only there are the statistical errors small enough to be interesting to a non-lattice 
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Figure 2: Rho mass vs. box size. 
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Figure 3: An Edinburgh plot for staggered fermions, from the MILC collaboration. 



audience. Most dynamical fermion simulations are unfortunately so noisy that it is hard to 
subject them to detailed questioning. 

Fig. U shows a plot of the rho mass as a function of the size of the simulation, for several 
values of the quark mass (or m w /m p ratio in the simulation) and lattice spacing (/? = 6.0 is 
a ~ 0.1 fm and (3 — 5.7 is about twice that).0 This picture shows that if the box has a diameter 
bigger than about 2 fm, the rho mass is little affected, but if the box is made smaller, the rho 
is "squeezed" and its mass rises. 

Next we look at an Edinburgh plot, Fig. The different plotting symbols correspond to 
different bare couplings or (equivalcntly) different lattice spacings. This plot shows large scaling 
violations: mass ratios from different lattice spacings do not lie on top of each other. We can 
expose the level of scaling violations by taking "sections" through the plot and plot tun /m p at 
fixed values of the quark mass (fixed m^/mp), vs. lattice spacing, in Fig. |j 
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Figure 4: "Sections" through the Edinburgh plot. 



Now for some examples of scaling tests in the chiral limit. (Extrapolating to the chiral limit 
is a whole can of worms on its own, but for now let's assume we can do it.) Fig. | shows the 
nucleon to rho mass ratio (at chiral limit) vs. lattice spacing (in units of 1/m p ) for staggeredQ 
and WilsonEl fermions. The "analytic" result is from strong coupling. The two curves are 
quadratic extrapolations to zero lattice spacing using different sets of points from the staggered 
data set. The burst is from a linear extrapolation to the Wilson data. The reason I show this 
figure is that one would like to know if the continuum limit of quenched spectroscopy "predicts" 
the real- world N/p mass ratio of 1.22 or not. The answer (unfortunately) depends on how the 
reader chooses to extrapolate. 

Another test0 is the ratio of the rho mass to the square root of the string tension, Fig. |^. 
Here the diamonds are staggered data and the crosses from the Wilson action. Scaling violations 
are large but the eye extrapolates to something close to data (the burst). 

Finally, despite Mr. Strange, very few authors have attempted to extrapolate to infinite 
volume, zero lattice spacing, and to physical quark masses, including the strange quark. One 
group which did, Butler et al.Jil produced Fig. [7| The squares are lattice data, the octagons are 
the real world. They look quite similar within errors. Unfortunately, to produce this picture, 
they had to build their own computer. 



3 Doing a Better Job— Maybe! 

The slow approach to scaling presents a practical problem for QCD simulations, since it means 
that one needs to work at small lattice spacing. This is expensive. The cost of a Monte Carlo 
simulation in a box of physical size L with lattice spacing a and quark mass m q scales roughly 
as 

(-n-r 2 (— r 3 m 

a a m a 
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Figure 5: Nucleoli to rho mass ratio (at chiral limit) vs. lattice spacing (in units of l/m p ). 
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Figure 6: Scaling test for the rho mass in terms of the string tension, with data points labeled 
as in Fig. 5. 
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Figure 7: Quenched approximation mass ratios from Ref. 13. 



where the 4 is just the number of sites, the 1-2 is the cost of "critical slowing down" -the extent 
to which successive configurations are correlated, and the 2-3 is the cost of inverting the fermion 
propagator, plus critical slowing down from the nearly massless pions. The problem is that one 
needs a big computer to do anything. 

However, all the simulations I described in the last lecture were done with a particular choice 
of lattice action: the plaquette gauge action, and either Wilson or staggered quarks. While 
those actions are the simplest ones to program, they are just particular arbitrary choices of bare 
actions. Can one invent a better lattice discretization, which has smaller scaling violations? 

People are trying many approaches. One could just write down a slightly more complicated 
action, include some parameters which can be tuned, do a spectroscopy calculation, and see if 
there is any improvement as the parameters are varied. The problem with this method is that 
it is like hunting for a needle in a multidimensional haystack- there are so many possible terms 
to add. One needs an organizing principle. 



3.1 Improvement based on naive dimensional analysis 

The simplest idea is to use the naive canonical dimensionality of operators to guide us in our 
choice of improvement. If we perform a naive Taylor expansion of a lattice operator like the 
plaquette, we find that it can be written as 

1 - ^Re TrU plaq = r TrF 2 u + a 2 [n £^ TtD^D^ + 

+0{a i ) (42) 

The expansion coefficients have a power series expansion in the coupling, rj = Aj + g 2 Bj + . . . 
and the expectation value of any operator T computed using the plaquette action will have an 
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Figure 8: The value of some parameter in a lattice action for which physical observables have 
no a n errors. The dotted line is the lowest order perturbative expectation. 



expansion 

(T(a)) = (T(0)) + O(a) + 0{g 2 a) + ... (43) 

Other loops have a similar expansion, with different coefficients. Now the idea is to take 
the lattice action to be a minimal subset of loops and systematically remove the a™ terms for 
physical observables order by order in n by taking the right linear combination of loops in the 
action. 

S = J2 c j°j ( 44 ) 

3 

with 

Cj = c° + g 2 c] + ... (45) 

This method was developed by Symanzik and co-workers&0 ten years ago. 

To visualize this technique, look at Fig. @. We imagine parameterizing the coefficients of 
various terms in the lattice action, which for a pure gauge theory could be a simple plaquette, 
a 1 x 2 closed loop, the square of the 1x2 loop, and so on, as some function of g 2 . "Tree- 
level improvement" involves specifying the value of the j-th coefficient Cj(g 2 ) at g 2 = 0. As we 
move away from g 2 = 0, the value of Cj(g 2 ) for which observables calculated using the lattice 
action have no errors through the specified exponent n (no a n errors) will trace out a trajectory 
in coupling constant space. For small g 2 , the variation should be describable by perturbation 
theory, Eq. [fB], but when g 2 gets large, we would not expect perturbation theory would be a 
good guide. 

The most commonly used "improved" fermion action is the "Sheikholeslami-Wohlert"lll or 
"clover" action, an order a 2 improved Wilson action. The original Wilson action has 0(a) errors 
in its vertices, Syy — S c + 0(a). This is corrected by making a field redefinition 

^{x)^^'{x)= + f 7„£VV> (46) 

$ix)^${x)= i>{x) + i i 1 ^D ll (47) 



Figure 9: The "clover term" . 



and the net result is an action with an extra lattice anomalous magnetic moment term, 



It is called the "clover" action because the lattice version of is the sum of paths shown in 



Studies performed at the time showed that this program did not improve scaling for the pure 
gauge theory (in the sense that the cost of simulating the more complicated action was greater 
than the savings from using a larger lattice spacing.) The whole program was re-awakened in 
the last few last years by Lepage and collaboratorsO and variations of this program give the 
most widely used "improved" lattice actions. 

3.2 Nonperturbative determination of coefficients 

Although I am breaking chronological order, the simplest approach to Symanzik improvement 
is the newest. The ideaEl is to force the lattice to obey various desirable identities to some order 
in a, by tuning parameters until the identities are satisfied by the simulations. That is, we try to 
find the solid line in Fig, || by doing simulations. Then use the action to calculate other things 
and test to see if scaling is improved. One example is the PCAC relation 



Ssw -7-^P{x)cr^F t _ lu ^(x) 



(48) 



Fig. § 



d„Al = 2m% + 0(a), 



(49) 



where the axial and pseudoscalar currents are just 




(50) 



and 



P a (x) 



(51) 
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Figure 10: Values of the quark mass as computed from the axial and pseudoscalar currents, 
using the Wilson action. The open and full symbols correspond to different boundary conditions 
on the gauge fields. 



(r a is an isospin index.) The PCAC relation for the quark mass is 

Now the idea is to take some Symanzik-improved action, with the improvement coefficients 
allowed to vary, and perform simulations in a little box with some particular choice of boundary 
conditions for the fields. Parameters which can be tuned include the csw i n the clover term 
i/Acswo-CfivFuv and ones used for more complicated expressions for the currents 

= Z A [(l + b A am q )A a ^ + c A ad ti P a (53) 

P a = Z P (l + b P am q )P a . (54) 

They are varied until the quark mass, defined in Eq. |5^, is independent of location in the box, 
or of the boundary conditions. Figs. [Io| and [ll] illustrate what can be done with this tuning 
procedure. It is still too soon for definitive tests of scaling with this procedure. 



3.3 Improving perturbation theory 

The older version of Symanzik improvement uses lattice perturbation theory to compute the 
coefficients of the operators in the action. The idea here is to find a new definition of g 2 for 
which the solid line in Fig. || is transformed into a straight line. (Compare Fig. |l2| .) 

Let's make a digression into lattice perturbation theory .0 It has three major uses. First, 
we need to relate lattice quantities (like matrix elements) to continuum ones: O cont (/i) = 
Z(fj,a, g(a))O latt (a). This happens because the renormalization of an operator is slightly dif- 
ferent in the two schemes. In perturbation theory Z has an expansion in powers of g 2 . Second, 
we can use perturbation theory to understand and check numerical calculations when the lattice 
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Figure 11: Same as previous figure, but now with improved action and operators. 
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Figure 12: Same as Fig. 8, but now with a redefined coupling constant g 2 R , to make the 
"improvement line" linear in g\. 



couplings are very small. Finally, one can use perturbative ideas to motivate nonperturbative 
improvement schemes.il 

Perturbation theory for lattice actions is just like any other kind of perturbation theory (only 
much messier). One expands the Lagrangian into a quadratic term and interaction terms and 
constructs the propagator from the quadratic terms: 

C= A fl {x)p^(x-y)A u {y) + gA 3 + ... (55) 
C + d. (56) 

For example, the gluon propagator in Feynman gauge for the Wilson action is 

ZV(<?) = ^ t: ^ , »• (57) 

E M (! - cos(g M o)) 

To do perturbation theory for any system (not just the lattice) one has to do three things: one 
has to fix the renormalization scheme (RS) (define a coupling), specify the scale at which the 
coupling is defined, and determine a numerical value for the coupling at that scale. All of these 
choices are arbitrary, and any perturbative calculation is intrinsically ambiguous. 
Any object which has a perturbative expansion can be written 

0(Q) = c + ci(Q//i, RS)a s (n, RS) + c 2 (Q/fx, RS)a s (n, RS) 2 + ... (58) 

In perturbative calculations we truncate the series after a fixed number of terms and implicitly 
assume that's all there is. The coefficients a(Q/fi, RS) and the coupling a s ([i,RS) depend on 
the renormalization scheme and choice of scale fi. The guiding rule of perturbation theoryB is 
"For a good choice of expansion the uncalculated higher order terms should be small." A bad 
choice has big coefficients. 

There are many ways to define a coupling: The most obvious is the bare coupling; as we will 
see shortly, it is a poor expansion parameter. Another possibility is to define the coupling from 
some physical observable. One popular choice is to use the heavy quark potential at very high 
momentum transfer to define 

V{ q ) = **C f *ML. (59) 
Q 

There are also several possibilities for picking a scale: One can use the bare coupling, then 
/i = l/a the lattice spacing. One can guess the scale or or play games just like in the continuum. 
One game is the Lepage-Mackenzie q* prescription: find the "typical" momentum transfer q* 
for a process involving a loop graph by pulling a s (q) out of the loop integral and set 

«.(?*) / A£(<?) - / A£(?H(g). (60) 

To find q* , write a s (q) = a s (p) + 61n(g 2 //i 2 )a s (/i) 2 + . . ., and similarly for a s (q*), insert these 
expressions into Eq. |6Fj and compare the a s (n) 2 terms, to get 

1%*) = J Aln(<z)£/ / A£. (61) 

This is the lattice analog of the Brodsky-Lepage-MackenzieEl prescription in continuum PT. 



Finally one must determine the coupling: If one uses the bare lattice coupling it is already 
known. Otherwise, one can compute it in terms of the bare coupling: 



a Ms( s / a ) = "o + (5.88 - 1.751ns)o$ + (43.41 - 21.891ns + 3.061n 2 s)ag + . . . (62) 

Or one can determine it from something one measures on the lattice, which has a perturbative 
expansion. For example 

1 47T 

- ln(~TrU plag ) = — a P (3.41/o)(l - 1.185a P ) (63) 
(to this order, ap — ay)- Does "improved perturbation theory" actually improve perturbative 



calculations? In many cases, yes: some examples are shown in Fig. 13 from B: On the upper 
left we see a calculation of the average link in Landau gauge, from simulations (octagons) and 
then from lowest-order perturbative calculations using the bare coupling (crosses) and ay and 
a TTs (diamonds and squares). In the upper right panel we see how the lattice prediction of 
an observable involving the 2 by 2 Wilson loop depends on the choice of momentum q* /a (at 
/3 = 6.2, a rather weak value of the bare coupling) in the running coupling constant. The burst 
is the value of the prescription of Eq. |6l]. In the lower panel are perturbative predictions the 
same observables as a function of lattice coupling. These pictures illustrate that perturbation 
theory in terms of the bare coupling does not work well, but that using other definitions for 
couplings, one can get much better agreement with the lattice "data" . 

Straight perturbative expansions by themselves for the commonly-used lattice actions are 
typically not very convergent. The culprit is the presence of U^s in the action. One might 
think that for weak coupling, one could expand 

$U%I> = + igaA + . . .]if> (64) 

and ignore the . . ., but the higher order term ip^g 2 a 2 A 2 ip generates the "tadpole graph" of 
Fig. [l4|. The UV divergence in the gluon loop ~ 1/a 2 cancels the a 2 in the vertex. The same 
thing happens at higher order, and the tadpoles add up to an effective a c n g 2n contribution. 
Parisill and later Lepage and Mackenzie^ suggested a heuristic way to deal with this problem: 
replace — ► Uq{1 + igaA) where uq, the "mean field term" or "tadpole improvement term" is 
introduced phenomenologically to sum the loops. Then one rewrites the action as 

^^"i^™ (65) 

where g 2 = gf at /uQ is the new expansion parameter. Is ufj = (Tr[/ p ; a q/3)? This choice is often 
used; it is by no means unique. 

A "standard action" (for this year, anyway) is the "tadpole-improved Liischer-Weisz0 ac- 
tion," composed of a 1 by 1, 1 by 2, and "twisted" loop (+x,+y,+z,-x,-y,-z), 
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PS = -/3[Tr(l x 1) - —5(1 + 0.48a s )Tr(l x 2) - 3 0.33a s TrU tw ] (66) 



with uo = (Tr[7pz ag /3) 1//4 and 3.068a s = — ln(TrJ7 p ; ag /3) determined self-consistently in the 
simulation. 
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Figure 13: Examples of "improved perturbation theory". 



Figure 14: The "tadpole diagram" . 



3.4 Fixed Point Actions 



Let's recall the question we were trying to answer in the previous sections: Can one find a 
trajectory in coupling constant space, along which the physics has no corrections to some desired 
order in a n or g m a n l Let's take the question one step further: Is there a trajectory in coupling 
constant space in which there are no corrections at all, for any n or ml 

To approach the answer, let's think about the connection between scaling and the properties 
of some arbitrary bare action, which we assume is defined with some UV cutoff a (which does 
not have to be a lattice cutoff, in principle). The action is characterized by an infinite number 
of coupling constants, {c}. (Many of them could be set to zero.) When the c's take on almost 
any arbitrary values, the typical scale for all physics will be the order of the cutoff: m ~ 1/a, 
correlation length £ ~ a. There will be strong cutoff effects. 

The best way to think about scaling is through the renormalization group.il Take the action 
with cutoff a and integrate out degrees of freedom to construct a new effective action with a new 
cutoff a' > a and fewer degrees of freedom. The physics at distance scales r > a is unaffected by 
the coarse-graining (assuming it is done exactly.) We can think of the effective actions as being 
similar to the original action, but with altered couplings. We can repeat this coarse-graining 
and generate new actions with new cutoffs. As we do, the coupling constants "flow:" 

S{a, cj) -> S(a', -> S(a", c'!) -»■ . . . (67) 

If under repeated blockings the system flows to a fixed point 

S(a n , c?) -> S(a n+1 ,c] +1 = c?) (68) 

then observables are independent of the cutoff a and in particular the correlation length £ must 
cither be zero or infinite. 

This can only happen if the original c's belong to a particular restricted set, called the 
"critical surface." It is easy to see that physics on the critical surface is universal: at long 
distances the effective theory is the action at the fixed point, to which all the couplings have 
flowed, regardless of their original bare values. In particular, physics at the fixed point is 
independent of the underlying lattice structure. 

But £ = oo is not £ large. Imagine tuning bare parameters close to the critical surface, 
but not on it. The system will flow towards the fixed point, then away from it. The flow 
lines in coupling constant space will asymptotically approach a particular trajectory, called the 
renormalized trajectory (RT), which connects (at £ = oo) with the fixed point. Along the 
renormalized trajectory, £ is finite. However, because it is connected to the fixed point, it shares 
the scaling properties of the fixed point-in particular, the complete absence of cutoff effects in 
the spectrum and in Green's functions. (To see this remarkable result, imagine doing QCD 
spectrum calculations with the original bare action with a cutoff equal to the Planck mass and 
then coarse graining. Now exchange the order of the two procedures. If this can be done without 
making any approximations the answer should be the same.) 

A Colorado analogy is useful for visualizing the critical surface and renormalized trajectory: 
think of the critical surface as the top of a high mountain ridge. The fixed point is a saddle 



RT 



/ 


\ 


z 1 


_\ . . . 


/ 


....A....... 







f = 5 



|=10 
CS 



Figure 15: A schematic picture of renormalization group flows along a one-dimensional critical 
surface, with the associated renormalized trajectory, and superimposed contours of constant 
correlation length. 

point on the ridge. A stone released on the ridge will roll to the saddle and come to rest. If it is 
not released exactly on the ridge, it will roll near to the saddle, then go down the gully leading 
away from it. For a cartoon, see Fig. [l5|. 

So the ultimate goal of "improvement programs" is to find a true perfect action, without 
cutoff effects, along the renormalized trajectory of some renormalization group transformation. 
At present, finding an RT has not been done in a convincing way for any renormalization group 
transformation. However, an action at the fixed point might also be an improved action, and 
fixed point actions really can be constructed and used. 

In lattice language, a bare action for QCD is described by one overall factor of j3 = 2N/g 2 
and arbitrary weights of various closed loops, 

2 AT 

5 3 

Asymptotic freedom is equivalent to the statement that the critical surface of any renormaliza- 
tion group transformation is at g 2 = 0. The location of a fixed point involves some relation 
among the Cj 's. 

A direct attack on the renormalized trajectory begins by finding a fixed point action. Imagine 
having a set of field variables {<fi} defined with a cutoff a. Introduce some coarse-grained variables 
{$} defined with respect to a new cutoff a', and integrate out the fine-grained variables to 
produce a new action 

e -/35'(*) = j #e -/3(T(#,,A)+S(0)) (70) 

where /3(T($, <ft) is the blocking kernel which functionally relates the coarse and fine variables. 
Integrating Eq. |7^is usually horribly complicated. However, P. Hasenfratz and F. NiedermayeiS 
noticed an amazing simplification for asymptotically free theories: Their critical surface is at 



= oo and in that limit Eq. [7(] becomes a steepest-descent relation 

S"($)=min((r($,^) + 5(^)) (71) 

which can be used to find the fixed point action 

S FP (^) = min((T($, 0) + S FP {4>))- (72) 

The program has been successfully carried out for d — 2 sigma models!! and for four-dimensional 
pure gauge theories.il These actions have two noteworthy properties: First, not only arc they 
classically perfect actions (they have no a n scaling violations for any n), but they are also 
one-loop quantum perfect: that is, as one moves out the renormalized trajectory, 

\s RT (g 2 ) = \(S FP + 0( ff 4 )). (73) 

g g 

Physically this happens because the original action has no irrelevant operators, and they are only 
generated through loop graphs. Thus these actions are an extreme realization of the Symanzik 
program. Second, because these actions are at the fixed point, they have scale invariant classical 
solutions. This fact can be used to define a topological charge operator on the lattice in a way 
which is consistent with the lattice action.S 

These actions are "engineered" in the following way: one picks a favorite blocking kernel, 
which has some free parameters, and solves Eq. ^2|, usually approximately at first. Then one 
tunes the parameters in the kernel to optimize the action for locality, and perhaps refines the 
solution. Now the action is used in simulations at finite correlation length (i.e. do simulations 
with a Boltzman factor exp(— (3S F p). Because of Eq. [73|, one believes that the FP action will 
be a good approximation to the perfect action on the RT; of course, only a numerical test can 
tell. As we will see in the next section, these actions perform very well. At this point in time no 
nonperturbative FP action which includes fermions has been tested, but most of the formalism 
is there. E§ 

3.5 Examples of "Improved" Spectroscopy 

1 would like to show some examples of the various versions of "improvement" , and remind you 
of the pictures at the end of the last chapter to contrast results from standard actions. 



Fig. 16 shows a plot of the string tension measured in systems of constant physical size 
(measured in units of 1/T C , the critical temperature for deconfinement), for SU(3) pure gauge 
theory. In the quenched approximation, with yfa ~ 440 MeV, T c = 275 MeV and 1/T C — 
0.7 fm. Simulations with the standard Wilson action are crosses, while the squares show FP 
action resultsEll and the octagons from the tadpole-improved Luscher-Weisz action.E3 The figure 
illustrates that it is hard to quantify improvement. There are no measurements with the Wilson 
action at small lattice spacing of of precisely the same observables that the "improvement people" 
measured. The best one can do is to take similar measurements (the diamonds) and attempt to 
compute the a = prediction for the observable we measured (the fancy cross at a = 0). This 
attempt lies on a straight line with the FP action data, hinting strongly that the FP action is 
indeed scaling. The FP action seems to have gained about a factor of three to four in lattice 
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Figure 16: The (square root of) the string tension in lattices of constant physical size L = 2/T c , 
but different lattice spacings (in units of 1/T C ). 



spacing, or a gain of (3 — 4) 6 compared to the plaquette action, according to Eq. [0], at a 
cost of a factor of 7 per site because it is more complicated to program. The tadpole-improved 
Liischer-Weisz action data lie lower than the FP action data and do not scale as well. As a — > 
the two actions should yield the same result; that is just universality at work. However, there 
is no guarantee that the approach to the continuum is monotonic. 

Fig. [It] shows the heavy quark-antiquark potential in SU(2) gauge theory, where V(r) and 
r are measured in the appropriate units of T c , the critical temperature for deconfincmcnt. The 
Wilson action is on the left and a FP action is on the right. The vertical displacements of the 
potentials are just there to separate them. Notice the large violations of rotational symmetry in 
the Wilson action data when the lattice spacing is a — 1/2T C which are considerably improved 
in the FP action results. 

Next we consider nonrelativistic QCD. A comparison of the quenched charmonium spectrum 
from Ref.0 using data from Ref.0 is shown in Fig. [l8| When the tadpole-improved L-W action 
is used to generate gauge configurations, the scaling window is pushed out to a ~ 0.4 fm for 
these observables. 

Now we turn to tests of quenched QCD for light quarks. The two actions which have been 
most extensively tested are the S-W action, with and without tadpole improvement, and an 
action called the D234(2/3) action, a higher-order variant of the S-W action.il Figs. Il9| and 
, are the analogs of Figs. 1 and § Diamond!! and plussesl are S-W actions, ordinary and 
tadpole- improved, squares are the D234(2/3) action. They appear to have about half the scaling 
violations as the standard actions but they don't remove all scaling violations. It's a bit hard to 
quantify the extent of improvement from these pictures because a chiral extrapolation is hidden 
in them. However, one can take one of the "sections" of Fig. ^ and overlay the new data on 
it, Fig. |l]. It looks like one can double the lattice spacing for an equivalent amount of scale 
violation. However, the extrapolation in a is not altogether clear. Fig. ^2] is the same data as 
Fig. Ell only plotted vs. a 2 , not a. All of the actions shown in these figures are supposed to 
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Figure 17: The heavy quark potential in SU(2) pure gauge theory measured in units of T c . (a) 
Wilson action (b) an FP action. 
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Figure 18: 5, P, and D states of charmonium computed on lattices with: a = 0.40 fm (improved 
action, (3 p i aq = 6.8); a — 0.33 fm (improved action, (3 p i aq = 7.1); a = 0.24 fm (improved action, 
fipiaq = 7.4); and a = 0.17 fm (Wilson action, fi = 5.7, from [43]), from Ref. 19. The dashed 
lines indicate the true masses. 
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Figure 19: Nucleon to rho mass ratio (at chiral limit) vs. lattice spacing (in units of l/m p ). 



have 0(a 2 ) (or better) scaling violations. Do the data look any straighter in Fig. E2 than in 



Fig. 21? 



3.6 The bottom line 

At the cost of enormous effort, one can do fairly high precision simulations of QCD in the 
quenched approximation with standard actions. The actions I have shown you appear to reduce 
the amount of computation required for pure gauge simulations from supercomputers to very 
large work stations, probably a gain of a few hundreds. All of the light quark data I showed 
actually came from supercomputers. According to Eq. ^l], a factor of 2 in the lattice spacing 
gains a factor of 64 in speed. The cost of either of the two improved actions I showed is about 
a factor of 8-10 times the fiducial staggered simulation. Improvement methods for fermions are 
a few years less mature than ones for pure gauge theory, and so the next time you hear a talk 
about the lattice, things will have changed for the better (maybe). 



4 SLAC Physics from the Lattice 

One of the major goals of lattice calculations is to provide hadronic matrix elements which either 
test QCD or can be used as inputs to test the standard model. In many cases the lattice numbers 
have uncertainties which are small enough that they are interesting to experimentalists. I want 
to give a survey of lattice calculations of matrix elements, and what better way at this summer 
school, than to recall science which was done here at SLAC, as the framework? 



4.1 Generic Matrix Element Calculations 

Most of the matrix elements measured on the lattice are expectation values of local operators 
composed of quark and gluon fields. The mechanical part of the lattice calculation begins 
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Figure 20: Rho mass scaling test with respect to the string tension. 
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Figure 21: mjv/m p vs am p at fixed quark mass (fixed 7n^/m p ). Interpolations of the S-W and 
D234(2/3) data were done by me. 
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Figure 22: m^ jm p vs (am p ) 2 at fixed quark mass (fixed m w /m p ). 



by writing down some Green's function which contains the local operator (call it J(x)) and 
somehow extracting the matrix element. For example, if one wanted (0\J(x)\h) one could look 
at the two-point function 

Cj O (t) = J2(°\ J (x,t)O(0,0)\0}. (74) 



Inserting a complete set of correctly normalized momentum eigenstates 

_ 1 ^ \A,p)(A,p\ 

L3 h 2Ea{p) 

and using translational invariance and going to large t gives 

___ mAt (0\J\A)(A\O\0) 



C JO (t) = e 



A second calculation of 

Coo(t) = J2(0\O(x,t)O(0,0)\0) = e 
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is needed to extract (0| J\A) by fitting two correlators with three parameters. 
Similarly, a matrix element (h\ J\h') can be gotten from 

Cab (t, t') = ^2(0\O A (t)J(x, t')O B (0)\0). 



(75) 



(76) 



(77) 



(78) 



by stretching the source and sink operators A and Ob far apart on the lattice, letting the lattice 
project out the lightest states, and then measuring and dividing out (0|Oa|^) and (0|Ob|/i). 

These lattice matrix elements are not yet the continuum matrix elements. The lattice is a 
UV regulator and changing from the lattice cutoff to a continuum regulator (like MS) introduces 
a shift 

(/lO^/i = l/a)|i) w = a D (l + ^-(Cj^ - C latt ) + . • .)(f\O latt (a)\i) + 0(a) + .... (79) 



The factor a D converts the dimensionless lattice number to its continuum result. The 0(a) 
corrections arise because the lattice operator might not be the continuum operator: df/dx = 
(f(x + a) — f(x))/a + 0(a). The C's are calculable in perturbation theory, and the "improved 
perturbation theory" described in the last section is often used to reduce the difference Cjj^ — 

Clatt- 



4.2 Structure Functions 

In the beginning there was deep inelastic scattering. The lattice knows about structure functions 
through their moments: 



/ dxx n ~ 2 F 2 (x, Q 2 ) = M n (Q 2 ) = V <Jv$P 
Jo 



(80) 



has a representation in terms of matrix elements of fairly complicated quark (for the nonsinglet 
structure function) or gluon bilinears: for quarks 
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OW-„ n = [^j Ti^ D m ~. D^q (81) 
£>, sp^Jp, s) ~ «(/)(p w ... Pf , n +...). (82) 

s 

D is a lattice covariant derivative, which is approximated by a finite difference. The Wilson 
coefficients cf^ are calculated in perturbation theory and depend on ji 2 /Q 2 as well as on the 



coupling constant g(/J,). The lattice calculation is done by sandwiching the operator in Eq. 78 

It is presently possible to calculate the two lowest moments of the proton structure function 
on the lattice. Two groupsEB presented results at this year's lattice conference. Fig. ^3] shows 
(x) from Ref.,@ a calculation done in quenched approximation. In this picture the massless 
quark limit is the left edge of the picture. There are several different lattice operators which 
serve as discretizations of the continuum operator, and the figure shows two possibilities. 

Unfortunately, the calculation is badly compromised by the quenched approximation. It 
shows (x) u — 0.38 and (x)d — 0.19, while in the real world we expect about 0.28 and 0.10 
respectively. In the computer there are no sea quarks, and their momentum is obviously picked 
up (at least partly) by the valence quarks. 

One nice feature about the lattice calculation is that the spin structure function can be 
calculated in essentially the same way; the operator O has an extra gamma-5 in it. RefS 
computed Am = 0.84 and Ad = —0.24 (in contrast to 0.92 and -0.34 in the real world). A plot 



vs quark mass is shown in Fig. 24 



There is no problem in principle which prevents extending these calculations to full QCD 
(with dynamical sea quarks). It will probably be very expensive to push beyond the lowest 
moments. 



4.3 Heavy Quark Physics 

Then there was the November revolution. Twenty years later, systems with one or more heavy 
quarks remain interesting objects for study. The lattice is no exception. Many groups study 



Figure 23: (x) for the proton (MOM scheme) from Ref. 35. The circles (boxes) correspond to 
different choices of lattice operators. The upper (lower) band of data represents the results for 
the up (down)-quark distribution. 
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Figure 24: Am (upper values) and Ad (lower values) for the proton from Ref. 35. 



GeV 
10.5 



10.0 



..05 $Q»«2S 



9.5- 



() 



.If. 



'So 3 Si 



— : Experiment 
□ : Fermilab (nr = 0) 
O: NRQCD (n/=0* 
• : NRQCD (n, = 2 
ffl : SCRl (n f = 2) 



O 



-2P 



ID 



-%#B1P 



'Pi 'D 2 



Figure 25: T Spectrum 



spectroscopy, decay amplitudes, form factors... with the goal of confronting both experiment 
and analytic theoretical models. 

There are several ways to study heavy quarks on the lattice. If the quark has infinite mass 
(the "static limit" ) its propagator is very simple: the quark is confined to one spatial location, 
and as it evolves in time, its color "twinkles." The propagator is just a product of link matrices 
going forward in time. 

One can simulate nonrelativistic quarks directly on the lattice.0 This has evolved into one 
of the most successful (and most elaborate) lattice programs. The idea is to write down lattice 
actions which are organized in an expansion of powers of the quark velocity and to systematically 
keep all the terms to some desired order. For example, one might write 

D 2 

S = ^[W t - — +p-B + ...]ip (83) 
2m 

including kinetic and magnetic moment terms, suitably (and artistically) discretized. Tadpole- 
improved perturbation theory is heavily used to set coefficients. Figs. ^H] and ^ show the 
Upsilon spectrum and its hyperfine splittings from various NRQCD calculations (from a recent 
summary by Shigemitsu0) . 

The main shortcoming of nonrelativistic QCD is of course that when the quark mass gets 
small the nonrelativistic approximation breaks down. For charmonium v/c ~ 0.3, so the method 
is less safe for this system than for the Upsilon. 

Finally, one can take relativistic lattice quarks and just make the mass heavy. If the quark 
mass gets too heavy (ma ~ 1) lattice artifacts dominate the calculation. For Wilson fcrmions, 
the dispersion relation breaks down: E(p) ~ mi +p 2 /2iri2 where m.2 7^ mi. The magnetic 
moment is governed by its own different mass, too. 

Another signal of difficulty is that all these formulations have their own pattern of scale 
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Figure 26: T Spin Splittings: Symbols have the same meaning as in Figure 25. 

violations. That is, nonrelativistic quarks and Wilson quarks approach their a — ► limits 
differently. This is often described in the literature by the statement that "the lattice spacing is 
different for different observables." For example, in one data set,@ the inverse lattice spacing (in 
MeV) is given as 2055 MeV from fitting the heavy qq potential, 2140 MeV from the rho mass, 
1800 MeV from the proton mass, and 2400 MeV from the Upsilon spectrum. These simulations 
are just sitting in the middle of figures like Fig. || with only one point, trying to guess where the 
left hand edge of the picture will be. This is a problem for calculations of B meson and baryon 
spectroscopy, where the heavy quarks might be treated nonrelativistically and the light quarks 
are relativistic. What observable should be used to set the overall scale?El 

One of the major uses of heavy quark systems by the lattice community is to try to calculate 
the strong coupling constant at Q 2 — M§. This topic deserves its own section. 



Now we are at the SLC and LEP. For some time now there have been claims that physics at the 
Z pole hints at a possible breakdown in the standard modelEl A key question in the discussion 
is whether or not the value of ctjjg inferred from the decay width of the Z is anomalously high 
relative to other determinations of the strong coupling (which are usually measured at lower Q 
and run to the Z pole). 

The most recent analysis of a s (Mz ) I am aware of is due to Erler and Langacker.0 Currently, 
^imeshape _ q j23(4) (2) (1) for the standard model Higgs mass range, where the first/ second/ 
third uncertainty is from inputs/ Higgs mass/ estimate of a 4 s terms. The central Higgs mass 
is assumed to be 300 GeV, and the second error is for M H = 1000 GeV (+), 60 GeV (-). For 
the SUSY Higgs mass range (60-150 GeV), one has the lower value aj^g = .121(4)(+1 - 0)(1). 
A global fit to all data gives 0.121(4)(1). Hinchcliffc in the same compilation quotes a global 
average of 0.118(3). 

The lattice can contribute to this question by predicting ajjg from low energy physics. The 
basic idea is simple: The lattice is a (peculiar) UV cutoff. A lattice mass fi = Ma plus an 
experimental mass M give a lattice spacing a = fi/M in fm. If one can measure some quantity 
related to a s at a scale Q ~ 1/a, one can then run the coupling constant out to the Z. 



4.4 a s (M z ) 
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Figure 27: Survey of aj^(M z ) from Ref. 42. 



The best (recent) lattice number, from Shigemitsu's Lattice 96 summary talkS 
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(84) 



where the first error includes both statistics and estimates of discretization errors, the second 
is due to uncertainties from the dynamical quark mass, and the third is from conversions of 
conventions. The lattice number is about one standard deviation below the pure Z-physics 
number. Lattice results are compared to other recent determinations of aj^g(Z) in Fig. |7[ a 
figure provided by P. Burrows. c3 

Two ways of calculating a s (Mz) from lattice have been proposed: The first is the "small loop 
method" .El This method uses the "improved perturbation theory" described in Chapter 2: One 
assumes that a version of perturbation theory can describe the behavior of short distance objects 
on the lattice: in particular, that the plaquette can be used to define ay(q = 3.41/a). With 
typical lattice spacings now in use, this gives the coupling at a momentum Qq = 8 — 10 GeV. 
One then converts the coupling to exj^g and runs out to the Z using the (published) three-loop 
beta function.0 

Usually, the lattice spacing is determined from the mass splittings of heavy QQ states. This 
is done because the mass differences between physical heavy quark states are nearly independent 
of the quark mass- for example, the S-P mass splitting of the -0 family is about 460 MeV, and it 
is about 440 MeV for the T. A second reason is that the mass splitting is believed to be much 
less sensitive to sea quark effects than light quark observables, and one can estimate the effects 
of sea quarks through simple potential models. The uncertainty in the lattice spacing is three 
to five per cent, but systematic effects are much greater (as we will see below). 

The coupling constant comes from Eq. |6^. The plaquette can be measured to exquisite 
accuracy (0.01 per cent is not atypical) and so the coupling constant is known essentially without 
error. However, the scale of the coupling is uncertain (due to the lattice spacing). 

The next problem is getting from lattice simulations, which are done with n/ — (quenched) 
or nt = 2 (but unphysical sea quark masses) to the real world of n/ = 3. Before simulations 



with dynamical fcrmions were available, the translation was done by running down in Q to a 
"typical gluonic scale" for the psi or the upsilon (a few hundred MeV) and then matching the 
coupling to the three- flavor coupling (in the spirit of effective field theories). This produced a 
rather low a s ~ 0.105. Now we have simulations at rif = 2 and can do better. Recall that in 
lowest order 



One measures l/a s in two simulations, one quenched, the other at n/ = 2, runs one measurement 
in Q to the Q of the other, then extrapolates 1/a linearly in nt to rif = 3. Then one can convert 
to M S and run away. 

Pictures like Fig. |2?j are not very useful when one wants to get a feel for the errors inherent 
in the lattice calculation. Instead, let's run our expectations for a s (Mz) down to the scale where 
the lattice simulations are done, and compare. Fig. ^8] does that. The squares are the results 
of simulations of charmed quarks and the octagons are from bottom quarks, both with nt = 0. 
The crosses and diamond are rif = 2 bottom and charm results. (The bursts show upsilon data 
when the 1S-2S mass difference gives a lattice spacing.) Note the horizontal error bars on the 
lattice data. Finally, the predicted rif = 3 coupling ap is shown in the fancy squares, with error 
bars now rotated because the convention is to quote an error in a s . The lower three lines in the 
picture (from top to bottom) are otjjg(Mz) = 0.118, 0.123, and 0.128 run down and converted 
to the lattice prescription. 

The two top lines are predictions for how quenched a should run. 

Now for the bad news. All of the rif = 2 data shown here were actually run on the same 
set of configurations. The bare couplings are the same, but the lattice spacings came out 
different. What is happening is that we are taking calculations at some lattice spacing and 
inferring a continuum numbers from them, but the lattice predictions have scale violations 
which are different. (The T calculations use nonrelativistic quarks, the ip calculations use heavy 
Wilson quarks.) Notice also that the bottom and charm quenched lattice spacings are different. 
This discrepancy is thought to be a failure of the quenched approximation: the characteristic 
momentum scale for binding in the ip and T are different, and because rif is not the real world 
value a runs incorrectly between the two scales. Said differently, in the quenched approximation, 
the spectrum of heavy quark bound states is different from the real world. 

There is a second method of determining a running coupling constant which actually allows 
one to see the running over a large range of scales. It goes by the name of the "Schrodinger 
functional," (referring to the fact that the authors study QCD in a little box with specified 
boundary conditions) but "coupling determined by varying the box size" would be a more 
descriptive title. It has been applied to quenched QCD but has not yet been extended to full 
QCT)jl and so it is not yet had an impact on phenomenology. This calculation does not use 
perturbation theory overtly. For a critical comparison of the two methods, se E. 



I have been hearing people talk about glueballs from psi decay for almost twenty years. TokiEl 
has summarized the experimental situation for glueballs. What do theorists expect for a spec- 
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4.5 Glueballs 
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Figure 28: Survey of ajjg(Q) at the scale where lattice simulations are actually done. 



trum? The problem is that any non-lattice model requires making uncontrolled approximations 
to get any kind of an answer: there are no obvious zeroth order solutions with small expansion 
parameters. The lattice is the only game in town for a first-principles calculation. 

People have been trying to measure the masses of the lightest glueballs (the scalar and the 
tensor) using lattice simulations for many years. The problem has proven to be very hard, for 
several reasons. 

Recall how we measure a mass from a correlation function (Eq. f!o| ). The problem with the 
scalar glueball is that the operator O has nonzero vacuum expectation value, and the correlation 
function approaches a constant at large t: 

lim C(t) -> \(0\O\p= 0)| 2 exp(-mi) + |(0|O|0)| 2 . (86) 



The statistical fluctuations on C (t) are given by Eq. |25| and we find after a short calculation 
that 

a -> —=. (87) 



Thus the signal to noise ratio collapses at large t like V N exp(-mt) due to the constant term. 

A partial cure for this problem is a good trial wave function O. While in principle the 
plaquette itself could be used, it is so dominated by ultraviolet fluctuations that it does not 
produce a good signal. Instead, people invent "fat links" which average the gauge field over 
several lattice spacings, and then make interpolating fields which are closed loops of these fat 
links. The lattice glueball is a smoke ring. 

The second problem is that lattice actions can have phase transitions at strong or inter- 
mediate coupling, which have nothing to do with the continuum limit, but mask continuum 
behavior El As an example of this, consider the gauge group SU(2), where a link matrix can be 
parameterized as U = lcos6* + ia ■ ns'mO, so Tr£7 = 2 cos 6*. Now consider a generalization of 
the Wilson action — S = (3TrU + 7(Tr£/) 2 . (this is a mixed fundamental-adjoint representation 
action). At 7 — > 00 Till — > ±1 and the gauge symmetry is broken down to Z(2). But Z(2) 



Figure 29: Phase transitions in the fundamental-adjoint plane. 



gauge theories have a first order phase transition. First order transitions are stable under per- 
turbations, and so the phase diagram of this theory, shown in Fig. has a line of first order 
transitions which terminate in a second order point. At the second order point some state with 
scalar quantum numbers becomes massless. However, now imagine that you are doing Monte 
Carlo along the 7 = line, that is, with the Wilson action. When you come near the critical 
point, any operator which couples to a scalar particle (like the one you are using to see the 
scalar glueball) will see the nearby transition and the lightest mass in the scalar channel will 
shrink. Once you are past the point of closest approach, the mass will rise again. Any scaling 
test which ignores the nearby singularity will lie to you. 

This scenario has been mapped out for SU(3), and the place of closest approach is at a 
Wilson coupling corresponding to a lattice spacing of 0.2 fm or so, meaning that very small 
lattice spacings are needed before one can extrapolate to zero lattice spacing. A summary of the 
situation is shown in Fig. |3^.E1 Here the quantity tq is the "Sommer radius" JE3 defined through 
the heavy quark force, F(r) = —dV(r)/dr, by r^F^ro) = —1.65. In the physical world of three 
colors and four flavors, ro = 0.5 fm. 

Finally, other arguments suggest that a small lattice spacing or a good approximation to an 
action on an RT are needed to for glueballs: the physical diameter of the glueball, as inferred 
from the size of the best interpolating field, is small, about 0.5 fm. Shafer and Shuryak0 have 
argued that the small size is due to instanton effects. Most lattice actions do bad things to 
instantons at large lattice spacing.il 

Two big simulations have carried calculations of the glueball mass close to the continuum 
limit: the UKQCD collaboration!! and a collaboration at IBM which built its own computer .0 
(The latter group is the one with the press release last December announcing the discovery of 
the glueball.) Their predictions in MeV are different and they each favor a different experi- 
mental candidate for the scalar glueball (the one which is closer to their prediction, of course). 
It is a useful object lesson because both groups say that their lattice numbers agree before 
extrapolation, but they extrapolate differently to a = 0. 
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Figure 30: Glueball mass vs ro with the Wilson action, from a summary picture in Ref. 49. 



The UKQCD group sees that the ratio to(0 ++ )/a/<t can be well fitted with a form b + ca 2 a (<r 
is the string tension) and a fit of this form to the lattice data of both groups gives m(0 ++ )/y/<r = 
3.64 ± 0.15. To turn this into MeV we need a in MeV units. One way is to take m p /y/(j and 
extrapolate that to a = using b + ca^/a. Averaging and putting 770 MeV for m p one gets 
y/a = 432 ±15 MeV, which is consistent with the usual estimate (from extracting the string 
tension from the heavy quark potential) of about 440 MeV. Using the total average they get 
m(0 ++ ) = 1572 ± 65 ± 55 MeV where the first error is statistical and the second comes from the 
scale. 

The IBM group, on the other hand, notices that m p a and m^a scale asymptotically, use 
the phi mass to predict quenched A-jg-g, then extrapolate m(0 ++ )/A = A + B(aA) 2 . They get 
1740(41) MeV from their data, when they analyze UKQCD data, they get 1625(94) MeV, and 
when they combine the data sets, they get 1707(64) MeV. 

A neutral reporter could get hurt here. It seems to me that the lattice prediction for the 
scalar glueball is 1600 ± 100 MeV, and that there are two experimental candidates for it, the 
/o(1500) and the /./(1710). 

Masses are not the end of the story. The IBM group has done two interesting recent calcu- 
lations related to glueballs, which strengthen their claim that the /j(1710) is the glueball. 

The first one of thcrnll was actually responsible for the press release. It is a calculation 
of the decay width of the glueball into pairs of pseudoscalars. This is done by computing an 
unamputated three point function on the lattice, with an assumed form for the vertex, whose 
magnitude is fitted. The result is shown in Fig. [3l|. The octagons are the results of the simulation 
and the diamonds show interpolations in the quark mass. The "experimental" points (squares) 
are from a partial wave analysis of isoscalar scalar resonances by Longacre and LindenbaumH 

The response of a member of the other side is that the slope of the straight line that one 
would put through the three experimental points is barely, if at all, compatible with the slope 
of the theoretical points. Since they argue theoretically for a straight line, the comparison of 
such slopes is a valid one. 
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Figure 31: Scalar glueball decay couplings from Ref. 54. 



If one of the experimental states is not a glueball, it is likely to be a 3 Po orbital excitation 
of quarks. Weingarten and LeeS are computing the mass of this state on the lattice and argue 
that it is lighter than 1700 MeV; in their picture the /o(1500) is an ss state. I have now said 
more than I know and will just refer you to recent discussions of the qucstion.0 

Both groups predict that the 2 ++ glueball is at about 2300 MeV. 

Can "improved actions" help the situation? Recently Peardon and MorningstarB imple- 
mented a clever method for beating the exponential signal-to-noise ratio: make the lattice 
spacing smaller in the time direction than in the space direction. Then the signal, which falls 
like exp(— matL t ) after L t lattice spacings, dies more slowly because at is reduced. Their picture 
of the glueball mass vs r<j. is shown in Fig. They are using the tadpole- improved Liischer- 
Weisz action. The pessimist notes the prominent dip in the middle of the curve; this action also 
has a lattice- artifact transition (somewhere); the optimist notes that the dip is much smaller 
than for the Wilson action and then the pessimist notes that there is no Wilson action data at 
large lattice spacing to compare. I think the jury is still out. 



4.6 The B Parameter of the B Meson 

And finally we are at BaBar. B — B mixing is parameterized by the ratio 
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= T b d^VQCDF{^)\VgV td \ 2 b(p){-(B\bl P ^ ~ K)dh P (l ~ l5)d\B)} (88) 
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Experiment is on the left; theory on the right. Moving into the long equation from the left, 
we see many known (more or less) parameters from phase space integrals or perturbative QCD 
calculations, then a combination of CKM matrix elements, followed by a four quark hadronic 
matrix elemcnt.il We would like to extract the CKM matrix element from the measurement 
of Xd (and its strange partner x s ). To do so we need to know the value of the object in the 
curly brackets, defined as 3/8Mbd and parameterized as ni 2 B fg d Bb d where Bf, d is the so-called 
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Figure 32: Glueball mass vs tq from Ref. 49, including large lattice spacing data. 



B-parameter, and fs is the B-meson decay constant 

(0\bjol5d\B) = f B m B . (89) 

Naive theory, which is expected to work well for the B system, suggests that Bb — 1 to good 
accuracy. Of course, the stakes are high and a good determination of Mm is needed to test the 
standard model. The lattice can do just that. 

In Eq. ^8] the coefficient which runs the effective interacion down from the W-boson 
scale to the QCD scale /i, and the matrix element M(/i) both depend on the QCD scale, and one 
often sees the renormalization group invariant quantities Mm = b(p)Albd(^) or Bm = 6(/i)i?bd(A*) 
quoted in the literature. 

Decay constants probe very simple properties of the wave function: in the nonrelativistic 
quark model 

fM = (90, 

where "0(0) is the qq wave function at the origin. For a heavy quark (Q) light quark (q) system 
■0(0) should become independent of the heavy quark's mass as the Q mass goes to infinity, and 
in that limit one can show in QCD that y / TOM~/ju approaches a constant. 

One way to compute the decay constant is to put a light quark and a heavy quark on the 
lattice and let them propagate. It is difficult to calculate Jb directly on present day lattices 
with relativistic lattice fermions because the lattice spacing is much greater than the b quark's 
Compton wavelength (or the UV cutoff is below m^). In this limit the b quark is strongly 
affected by lattice artifacts as it propagates. However, one can make infinite on the lattice 
and determine the combination y/msfB in the limit. Then one can extrapolate down to the 
B mass and see if the two extrapolations up and down give the same result. (Nonrelativistic b 
quarks can solve this problem in principle, but the problem of setting the lattice spacing between 
light and nonrelativistic quarks has prevented workers from quoting a useful decay constant from 
these simulations.) 
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Figure 33: Pseudoscalar mesor^decay constant vs 1/M, from Ref. 59 



Among the many lattice decay constant calculations, the one of Ref! stands out in my mind 
for being the most complete. These authors did careful quenched simulations at many values of 
the lattice spacing, which allows one to extrapolate to the continuum limit by brute force. They 
have also done a less complete set of simulations which include light dynamical quarks, which 
should give some idea of the accuracy of the quenched approximation. 

The analysis of all this data is quite involved. One begins with a set of lattice decay constants 
measured in lattice units, from simulations done with heavy quarks which are probably too light 
and light quarks which are certainly too heavy. One has to interpolate or extrapolate the heavy 
quark masses to their real world values, extrapolate the light quarks down in mass to their 
physical values, and finally try to extrapolate to a — > 0. It is not always obvious how to do 
this. Complicating everything are the lattice artifacts in the fermion and gauge actions, and the 
lattice-to-continuum renormalization factors as in Eq. [7^. 

The (still preliminary) results of Ref.El are shown in Figs. |33] and 34. The Nf = 2 dynamical 
fermion data in Fig. |^ have moved around a bit in the past year and may not have settled 
down yet. 

The numerical results of Ref.il are: 



f B = 166(11)(28)(14) f D = 196(9)(14)(8) 
f B , = 181(10)(36)(18) f Ds = 211(7)(25)(11) 

^ = 1.10(2)(5)(8) £t = 1.09(2)(5)(5) 

}B JD 

where the first error includes statistical errors and systematic effects of changing fitting ranges; 
the second, other errors within the quenched approximation; the third, an estimate of the 
quenching error. Decay constants are in MeV. 

Note that the error bars for the B system are small enough to be phenomenologically inter- 
esting. The Particle Data Group determination of CKM matrix elements, which does not 
include this data, says, "Using B B J Bi = (1.2±0.2)(173±40 MeV) 2 ..., \V* b V td \ 2 = 0.009±0.003, 
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Figure 34: fs vs. a from Ref. 59. Octagons are quenched data; crosses, Np = 2. The solid line 
is a linear fit to all quenched points; the dashed line is a a constant fit to the three quenched 
points with a < 0.5 GeV -1 . The extrapolated values at a = are indicated by bursts. The 
scale is set by /„• = 132 MeV throughout. 

where the error bar comes primarily from the theoretical uncertainty in the hadronic matrix el- 
ements." 

Can we trust these numbers? Lattice calculations have been predicting fu s ~ 200 MeV 
for about eight years. The central values have changed very little, while the uncertainties have 
decreased. So far four experiments have reported measurements of this quantity. The most 
recent is Fermilab E653 CollaboratiorH f Ds = 194(35_K20)(14) MeV. The older numbers, with 
bigger errors, were 238(47) (21) (43) from WA75 (1993)0; 344(37) (52) (42) from CLEO (1994)0; 
430 (+150 -130) (40) from BES (1995)0 

Now back to the mixing problem. On the lattice, one could measure the decay constants 
and B parameter separately and combine them after extrapolation, or measure M directly and 
extrapolate it. In principle the numbers should be the same, but in practice they will not be. 

A recent calculation illustrates this point.i The authors computed the four fermion oper- 
ators Mb s and Mtd directly on the lattice. Fig. |3^ shows the behavior of M as a function of 
hadron mass at one of their lattice spacings. 

The ratio r sc i = M^s/Mm is presumably much less sensitive to lattice spacing or to quark 
mass extrapolation. The authors' result for the lattice spacing dependence of this ratio is shown 
in Fig. |36|, along with an extrapolation to zero lattice spacing. They find 1.54 ± .13 ± .32 from 
their direct method, compared to r sc i « 1.32 ± .23 from separate extrapolations of the decay 
constants and the B parameter. (They measure equal B parameters for strange and nonstrange 
B mesons, B{fi) = 1.02(13) for fi = 2 GeV.) 

It looks like SU(3) breaking is fairly large, and if that is so, it looks like the parameter x s , 
the strange analog of Eq. might be about 20, unmeasurably largeS 
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Figure 35: Mh c (octagons) and Mh. s (squares) as a function of the inverse heavy-down(strangc) 
meson mass, at (3 = 6.0. The dashed line shows the effect of the lightest points on the fit. 
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Figure 36: The SU(3) flavor breaking ratio Mb s /Mbd versus the lattice spacing a. The points 
denoted by crosses were used in the fit (solid line). The burst shows the extrapolation to a = 0. 



5 Conclusions 



Lattice methods have arrived. There are so many lattice calculations of different matrix elements 
that it is impossible to describe them all, and in many cases the quality of the results is very 
high. One can see plots showing extrapolations in lattice spacing which show that the control of 
lattice spacing has become good enough to make continuum predictions with small uncertainties. 
Calculations with dynamical fermions and a small lattice spacing are still nearly impossibly 
expensive to perform, and "quenching" remains the dominant unknown in all lattice matrix 
element calculations. 

There are two major tasks facing lattice experts. I believe that all the people in our field 
would agree that the first problem is to reduce the computational burden, so that we can do 
more realistic simulations with smaller computer resources. I have illustrated several of the 
approaches people are using to attack this problem. I believe that some of them have been 
shown to be successful, and that "improvement" will continue to improve. 

There is a second question for lattice people, which I have not discussed, but I will mention 
at the end: Is there a continuum phenomenology of light hadron structure or confinement, which 
can be justified from lattice simulations? The motivation for asking this question is that there are 
many processes which cannot be easily addressed via the lattice, but for which a QCD prediction 
ought to exist. For examples of such questions, see the talk of Bjorken in this conference ,0 or 
Shuryak's article.Q Few lattice people are thinking about this question. Part of the lattice 
community spends its time looking for "structure" in Monte Carlo-generated configurations of 
gauge fields: instantons, monopoles, . . .. This effort is not part of the mainstream because the 
techniques either involve gauge fixing (and so it is not clear whether what is being seen is just 
an artifact of a particular gauge), or they involve arbitrary decisions during the search (perform 
a certain number of processing operations, no more, no less). To answer this question requires 
new ideas and a controlled approach to simulation data. Will an answer be found? 
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